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MEDIA H ANDLING SYSTEM 
This invention relates to media handling systems. 

Media handling systems such as so-called non-linear editing systems are widely used 

5 in the production of media (e.g. audio and video) material for broadcast and other 
applications. These systems are generally based on a computer workstation such as a general 
purpose personal computer, and allow sections of material, such as video and audio "clips", 
to be edited together into an output media product. Unlike previous (linear) editing systems, 
the clips can be handled, trimmed and assembled into the output media product in any order; 

10 that is to say, the editing operations themselves are not tied to the temporal order in which 
the material was originally captured. 

In such an editing system, a computer display will often show screen icons or other 
representations of the available media clips, so that the operator can select the ones that he 
wishes to include in the output media product. A problem, however, is that a production 

15 may involve very large numbers of clips, so that it is easy for the operator to lose track of 
which clip is which. Or the operator may wish to bring in a clip from an archive store, but it 
may take some considerable time to locate the correct one. 

This problem is alleviated to some extent by the use of so-called "picture stamps" 
which are small on-screen representations, associated with each clip icon, of what is 

20 perceived to be a representative picture from that clip. For example, the first, last or middle 
picture from the clip may be used as the picture stamp. 

It is a constant aim in this technical field to increase the efficiency of use of media 
handling systems. 

This invention provides a media handling system in which candidate video sequences 
25 are displayed on a display screen in schematic form for selection by a user, the system 
comprising: means for detecting human faces in the candidate video sequences; a display 
screen for displaying representations of the candidate video sequences for selection by a 
user, each representation including one or more images representing human faces derived 
from the respective video sequences; and a user control for defining a set of one or more of 
30 the video sequences. 

The invention provides an elegant technique for improving the efficiency of use of a 
media handling system such as an editing system, in that each candidate video sequence 
containing a human face has a representation of at least one face displayed with it, for 
example as a representative picture stamp, on the display screen. This means that the 
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operator can more easily recognise the content of the candidate video sequence from the face 
provided in the representative picture stamp. 

The invention is particularly suited for use in editing systems, in which the set of one 
or more of the video sequences is an ordered edited set forming an output media product. 
5 The invention also provides a method of media handling in which candidate video 

sequences are displayed on a display screen in schematic form for selection by a user, the 
method comprising the steps of: detecting human faces in the candidate video sequences; 
displaying on a display screen representations of the candidate video sequences for selection 
by a user, each representation including one or more images representing human faces 
10 derived from the respective video sequences; and providing a user control for defining a set 
of one or more of the video sequences. 

Further respective aspects and features of the invention are defined in the appended 

claims. 

Embodiments of the invention will now be described, by way of example only, with 
15 reference to the accompanying drawings, throughout which like parts are defined by like 
numerals, and in which: 

Figure 1 is a schematic diagram of a general purpose computer system for use as a 
face detection system and/or a non-linear editing system; 

Figure 2 is a schematic diagram of a video camera-recorder (camcorder) using face 
20 detection; 

Figure 3 is a schematic diagram illustrating a training process; 
Figure 4 is a schematic diagram illustrating a detection process; 
Figure 5 schematically illustrates a feature histogram; 

Figure 6 schematically illustrates a sampling process to generate eigenblocks; 
25 Figures 7 and 8 schematically illustrates sets of eigenblocks; 

Figure 9 schematically illustrates a process to build a histogram representing a block 
position; 

Figure 10 schematically illustrates the generation of a histogram bin number; 
Figure 1 1 schematically illustrates the calculation of a face probability; 
30 Figures 12a to 12f are schematic examples of histograms generated using the above 

methods; 

Figures 13a to 13g schematically illustrate so-called multiscale face detection; 
Figure 14 schematically illustrates a face tracking algorithm; 
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Figures 15a and 15b schematically illustrate the derivation of a search area used for 

skin colour detection; 

Figure 16 schematically illustrates a mask applied to skin colour detection; 
Figures 17a to 17c schematically illustrate the use of the mask of Figure 16; 
Figure 18 is a schematic distance map; 

Figures 19a to 19c schematically illustrate the use of face tracking when applied to a 
video scene; 

Figure 20 schematically illustrates a display screen of anon-linear editing system; 
Figures 21a and 21b schematically illustrate clip icons; Figures 22a to 22c 
schematically illustrate a gradient pre-processing technique; 

Figure 23 schematically illustrates a video conferencing system; 

Figures 24 and 25 schematically illustrate a video conferencing system in greater 

detail; 

Figure 26 is a flowchart schematically illustrating one mode of operation of the 

system of Figures 23 to 25; 

Figures 27a and 27b are example images relating to the flowchart of Figure 26; 
Figure 28 is a flowchart schematically illustrating another mode of operation of the 

system of Figures 23 to 25; 

Figures 29 and 30 are example images relating to the flowchart of Figure 28; 

Figure 31 is a flowchart schematically illustrating another mode of operation of the 

system of Figures 23 to 25; 

Figure 32 is an example image relating to the flowchart of Figure 31; and 

Figures 33 and 34 are flowcharts schematically illustrating further modes of 

operation of the system of Figures 23 to 25; 

Figure 1 is a schematic diagram of a general purpose computer system for use as a 
face detection system and/or a non-linear editing system. The computer system comprises a 
processing until 10 having (amongst other conventional components) a central processing 
unit (CPU) 20, memory such as a random access memory (RAM) 30 and non-volatile 
storage such as a disc drive 40. The computer system may be connected to a network 50 
, such as a local area network or the Internet (or both). A keyboard 60, mouse or other user 
input device 70 and display screen 80 are also provided. The skilled man will appreciate that 
a general purpose computer system may include many other conventional parts which need 
not be described here. 
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Figure 2 is a schematic diagram of a video camera-recorder (camcorder) using face 
detection. The camcorder 100 comprises a lens 110 which focuses an image onto a charge 
coupled device (CCD) image capture device 120. The resulting image in electronic form is 
processed by image processing logic 130 for recording on a recording medium such as a tape 
cassette 140. The images captured by the device 120 are also displayed on a user display 
150 which may be viewed through an eyepiece 160. 

To capture sounds associated with the images, one or more microphones are used. 
These may be external microphones, in the sense that they are connected to the camcorder by 
a flexible cable, or maybe mounted on the camcorder body itself. Analogue audio signals 
from the microphone (s) are processed by an audio processing arrangement 170 to produce 
appropriate audio signals for recording on the storage medium 140. 

It is noted that the video and audio signals may be recorded on the storage medium 
140 in either digital form or analogue form, or even in both forms. Thus, the image 
processing arrangement 130 and the audio processing arrangement 170 may include a stage 
of analogue to digital conversion. 

The camcorder user is able to control aspects of the lens 110's performance by user 
controls 180 which influence a lens control arrangement 190 to send electrical control 
signals 200 to the lens 110. Typically, attributes such as focus and zoom are controlled in 
this way, but the lens aperture or other attributes may also be controlled by the user. 

Two fiirther user controls are schematically illustrated. A push button 210 is 
provided to initiate and stop recording onto the recording medium 140. For example, one 
push of the control 210 may start recording and another push may stop recording, or the 
control may need to be held in a pushed state for recording to take place, or one push may 
start recording for a certain timed period, for example five seconds. In any of these 
arrangements, it is technologically very straightforward to establish from the camcorder's 
record operation where the beginning and end of each "shot" (continuous period of 
recording) occurs. 

The other user control shown schematically in Figure 2 is a "good shot marker" 
(GSM) 220, which may be operated by the user to cause "metadata" (associated data) to be 
stored in connection with the video and audio material on the recording medium 140, 
indicating that this particular shot was subjectively considered by the operator to be "good" 
in some respect (for example, the actors performed particularly well; the news reporter 
pronounced each word correctly; and so on). 
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The metadata may be recorded in some spare capacity (e.g. "user data") on the 
recording medium 140, depending on the particular format and standard in use. 
Alternatively, the metadata can be stored on a separate storage medium such as a removable 
MemoryStick R ™ memory (not shown), or the metadata could be stored on an external 
5 database (not shown), for example being communicated to such a database by a wireless link 
(not shown). The metadata can include not only the GSM information but also shot 
boundaries, lens attributes, alphanumeric information input by a user (e.g. on a keyboard - 
not shown), geographical position information from a global positioning system receiver (not 
shown) and so on. 

10 So far, the description has covered a metadata-enabled camcorder. Now, the way in 

which face detection may be applied to such a camcorder will be described. 

The camcorder includes a face detector arrangement 230. Appropriate arrangements 
will be described in much greater detail below, but for this part of the description it is 
sufficient to say that the face detector arrangement 230 receives images from the image 

15 processing arrangement 130 and detects, or attempts to detect, whether such images contain 
one or more faces. The face detector may output face detection data which could be in the 
form of a "yes/no" flag or maybe more detailed in that the data could include the image co- 
ordinates of the faces, such as the co-ordinates of eye positions within each detected face. 
This information may be treated as another type of metadata and stored in any of the other 

20 formats described above. 

As described below, face detection may be assisted by using other types of metadata 
within the detection process. For example, the face detector 230 receives a control signal 
from the lens control arrangement 190 to indicate the current focus and zoom settings of the 
lens 110. These can assist the face detector by giving an initial indication of the expected 

25 image size of any faces that may be present in the foreground of the image. In this regard, it 
is noted that the focus and zoom settings between them define the expected separation 
between the camcorder 100 and a person being filmed, and also the magnification of the lens 
110. From these two attributes, based upon an average face size, it is possible to calculate 
the expected size (in pixels) of a face in the resulting image data. 

30 A conventional (known) speech detector 240 receives audio information from the 

audio processing arrangement 170 and detects the presence of speech in such audio 
information. The presence of speech may be an indicator that the likelihood of a face being 
present in the corresponding images is higher than if no speech is detected. 



WO 2004/051656 



6 



PCT/GB2003/005188 



Finally, the GSM information 220 and shot information (from the control 210) are 
supplied to the face detector 230, to indicate shot boundaries and those shots considered to 
be most usefid by the user. 

Of course, if the camcorder is based upon the analogue recording technique, further 
5 analogue to digital converters (ADCs) may be required to handle the image and audio 
information. 

The present embodiment uses a face detection technique arranged as two phases. 
Figure 3 is a schematic diagram illustrating a training phase, and Figure 4 is a schematic 
diagram illustrating a detection phase. 

10 Unlike some previously proposed face detection methods (see References 4 and 5 

below), the present method is based on modelling the face in parts instead of as a whole. The 
parts can either be blocks centred over the assumed positions of the facial features (so-called 
"selective sampling") or blocks sampled at regular intervals over the face (so-called "regular 
sampling"). The present description will cover primarily regular sampling, as this was found 

15 in empirical tests to give the better results. 

In the training phase, an analysis process is applied to a set of images known to 
contain faces, and (optionally) another set of images ("nonface images") known not to 
contain faces. The analysis process builds a mathematical model of facial and nonfacial 
features, against which a test image can later be compared (in the detection phase). 

20 So, to build the mathematical model (the training process 310 of Figure 3), the basic 

steps are as follows: 

1. From a set 300 of face images normalised to have the same eye positions, each face is 
sampled regularly into small blocks. 

2. Attributes are calculated for each block; these attributes are explained further below. 
25 3. The attributes are quantised to a manageable number of different values. 

4. The quantised attributes are then combined to generate a single quantised value in 
respect of that block position. 

5. The single quantised value is then recorded as an entry in a histogram, such as the 
schematic histogram of Figure 5. The collective histogram information 320 in respect of all 

30 of the block positions in all of the training images forms the foundation of the mathematical 
model of the facial features. 

One such histogram is prepared for each possible block position, by repeating the 
above steps in respect of a large number of test face images. The test data are described 
further in Appendix A below. So, in a system which uses an array of 8 x 8 blocks, 64 
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histograms are prepared. In a later part of the processing, a test quantised attribute is 
compared with the histogram data; the fact thai a whole histogram is used to model the data 
means that no assumptions have to be made about whether it follows a parameterised 
distribution, e.g. Gaussian or otherwise. To save data storage space (if needed), histograms 
5 which are similar can be merged so that the same histogram can be reused for different block 
positions. 

In the detection phase, to apply the face detector to a test image 350, successive 
windows in the test image are processed 340 as follows: 

6. The window is sampled regularly as a series of blocks, and attributes in respect of 
1 o each block are calculated and quantised as in stages 1 -4 above. 

7. Corresponding "probabilities" for the quantised attribute values for each block 
position are looked up from the corresponding histograms. That is to say, for each block 
position, a respective quantised attribute is generated and is compared with a histogram 
previously generated in respect of that block position. The way in which the histograms give 

1 5 rise to "probability" data will be described below. 

8. All the probabilities obtained above are multiplied together to form a final probability 
which is compared against a threshold in order to classify the window as "face" or 
"nonface". It will be appreciated that the detection result of "face" or "nonface" is a 
probability-based measure rather than an absolute detection. Sometimes, an image not 

20 containing a face may be wrongly detected as "face", a so-called false positive. At other 
times, an image containing a face may be wrongly detected as "nonface", a so-called false 
negative. It is an aim of any face detection system to reduce the proportion of false positives 
and the proportion of false negatives, but it is of course understood that to reduce these 
proportions to zero is difficult, if not impossible, with current technology. 

As mentioned above, in the framing phase, a set of "nonface" images can be used to 
generate a corresponding set of "nonface" histograms. Then, to achieve detection of a face, 
the "probability" produced from the nonface histograms may be compared with a separate 
threshold, so that the probability has to be under the threshold for the test window to contain 
a face. Alternatively, the ratio of the face probability to the nonface probability could be 

30 compared with a threshold. 

Extra training data may be generated by applying "synthetic variations" 330 to the 
original ttaining set, such as variations in position, orientation, size, aspect ratio, background 
scenery, lighting intensity and frequency content. 



25 
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The derivation of attributes and their quantisation will now be described. In the 
present technique, attributes are measured with respect to so-called eigenblocks, which are 
core blocks (or eigenvectors) representing different types of block which may be present in 
the windowed image. The generation of eigenblocks will first be described with reference to 
5 Figure 6. 

Eigenblock creation 

The attributes in the present embodiment are based on so-called eigenblocks. The 
eigenblocks were designed to have good representational ability of the blocks in the training 
10 set. Therefore, they were created by performing principal component analysis on a large set 
of blocks from the training set. This process is shown schematically in Figure 6 and 
described in more detail in Appendix B. 

Training the System 

15 Experiments were performed with two different sets of training blocks. 

Eigenblock set I 

Initially, a set of blocks were used that were taken from 25 face images in the training 
set. The 16x16 blocks were sampled every 16 pixels and so were non-overlapping. This 
20 sampling is shown in Figure 6. As can be seen, 16 blocks are generated from each 64x64 
training image. This leads to a total of 400 training blocks overall. 

The first 10 eigenblocks generated from these training blocks are shown in Figure 7. 

Eigenblock set II 

25 A second set of eigenblocks was generated from a much larger set of training blocks. 

These blocks were taken from 500 face images in the training set In this case, the 16x16 
blocks were sampled every 8 pixels and so overlapped by 8 pixels. This generated 49 blocks 
from each 64x64 training image and led to a total of 24,500 training blocks. 

The first 12 eigenblocks generated from these training blocks are shown in Figure 8. 

30 Empirical results show that eigenblock set II gives slightly better results than set I. 

This is because it is calculated from a larger set of training blocks taken from face images, 
and so is perceived to be better at representing the variations in faces. However, the 
improvement in performance is not large. 
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Building the Histograms 

A histogram was built for each sampled block position within the 64x64 face image. 
The number of histograms depends on the block spacing. For example, for block spacing of 
16 pixels, there are 16 possible block positions and thus 16 histograms are used. 
5 The process used to build a histogram representing a single block position is shown 

in Figure 9. The histograms are created using a large training set 400 of M face images. For 
each face image, the process comprises: 

• Extracting 41 0 the relevant block from a position (i j) in the face image. 

• Calculating the eigenblock-based attributes for the block, and determining the relevant 
10 bin number 420 from these attributes. 

• Incrementing the relevant bin number in the histogram 430. 

This process is repeated for each of M images in the training set, to create a 
histogram that gives a good representation of the distribution of frequency of occurrence of 
the attributes. Ideally, M is very large, e.g. several thousand. This can more easily be 
15 achieved by using a training set made up of a set of original faces and several hundred 
synthetic variations of each original face. 

Generating the histogram bin number 

A histogram bin number is generated from a given block using the following process, 
20 as shown in Figure 10. The 16x16 block 440 is extracted from the 64x64 window or face 
image. The block is projected onto the set 450 of A eigenblocks to generate a set of 
"eigenblock weights". These eigenblock weights are the "attributes" used in this 
implementation. They have a range of -1 to +1. This process is described in more detail in 
Appendix B. Each weight is quantised into a fixed number of levels, L, to produce a set of 
25 quantised attributes 470, w„i = l.A. The quantised weights are combined into a single 

value as follows: 

h = w l L A ~ X + w 2 L A ~ 2 + w 3 L A ~ 3 + ... + w^L 1 + w A L° 
where the value generated, h, is the histogram bin number 480. Note that the total number of 

bins in the histogram is given by L A . 
30 The bin "contents", i.e. the frequency of occurrence of the set of attributes giving rise 

to that bin number, may be considered to be a probability value if it is divided by the number 
of training images M. However, because the probabilities are compared with a threshold, 
there is in fact no need to divide through by M as this value would cancel out in the 
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calculations. So, in the following discussions, the bin "contents" will be referred to as 
"probability values", and treated as though they are probability values, even though in a strict 
sense they are in fact frequencies of occurrence. 

The above process is used both in the training phase and in the detection phase. 

5 

Face Detection Phase 

The face detection process involves sampling the test image with a moving 64x64 
window and calculating a face probability at each window position. 

The calculation of the face probability is shown in Figure 1 1 . For each block position 
10 in the window, the block's bin number 490 is calculated as described in the previous section. 
Using the appropriate histogram 500 for the position of the block, each bin number is looked 
up and the probability 510 of that bin number is determined. The sum 520 of the logs of 
these probabilities is then calculated across all the blocks to generate a face probability 
value, P faee (otherwise referred to as a log likelihood value). 

15 This process generates a probability "map" for the entire test image. In other words, a 

probability value is derived in respect of each possible window centre position across the 
image. The combination of all of these probability values into a rectangular (or whatever) 
shaped array is then considered to be a probability "map" corresponding to that image. 

This map is then inverted, so that the process of finding a face involves finding 

20 minima in the inverted map. A so-called distance-based technique is used. This technique 
can be summarised as follows: The map (pixel) position with the smallest value in the 
inverted probability map is chosen. If this value is larger than a threshold (TD), no more 
faces are chosen. This is the termination criterion. Otherwise a face-sized block 
corresponding to the chosen centre pixel position is blanked out (i.e. omitted from the 

25 following calculations) and the candidate face position finding procedure is repeated on the 
rest of the image until the termination criterion is reached. 

Nonface method 

The nonface model comprises an additional set of histograms which represent the 
30 probability distribution of attributes in nonface images. The histograms are created in exactly 
the same way as for the face model, except that the training images contain examples of 
nonfaces instead of faces. 
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During detection, two log probability values are computed, one using the face model 
and one using the nonface model. These are then combined by simply subtracting the 
nonface probability from the face probability: 
p - P — P 

r combined ~ * face 1 nonface 

P combined is then used instead of P faee to produce the probability map (before 
inversion). 

Note that the reason that P nonfoce is subtracted from P faee is because these are log 
probability values. 



10 Histogram Examples 

Figures 12a to 12f show some examples of histograms generated by the training 

process described above. 

Figures 12a, 12b and 12c are derived from a training set of face images, and Figures 
12d, 12e and 12f are derived from a training set of nonface images. In particular: 

15 





Face histograms 


Nonface histograms 


Whole histogram 


Figure 12a 


Figure 12d 


Zoomed onto the main peaks at about h=l 500 


Figure 12b 


Figure 12e 


A further zoom onto the region about h=1570 


Figure 12c 


Figure 12f 



It can clearly be seen that the peaks are in different places in the face histogram and 
the nonface histograms. 

20 

Multiscale face detection 

In order to detect faces of different sizes in the test image, the test image is scaled by 
a range of factors and a distance (i.e. probability) map is produced for each scale. In Figures 
13a to 13c the images and their corresponding distance maps are shown at three different 
25 scales. The method gives the best response (highest probability, or niinimum distance) for 
the large (central) subject at the smallest scale (Fig 13a) and better responses for the smaller 
subject (to the left of the main figure) at the larger scales. (A darker colour on the map 
represents a lower value in the inverted map, or in other words a higher probability of there 
being a face).Candidate face positions are extracted across different scales by first finding 
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the position which gives the best response over all scales. That is to say, the highest 
probability (lowest distance) is established amongst all of the probability maps at all of the 
scales. This candidate position is the first to be labelled as a face. The window centred over 
that face position is then blanked out from the probability map at each scale. The size of the 

5 window blanked out is proportional to the scale of the probability map. 

Examples of this scaled blanking-out process are shown in Figures 13a to 13c. In 
particular, the highest probability across all the maps is found at the left hand side of the 
largest scale map (Figure 13c). An area 530 corresponding to the presumed size of a face is 
blanked off in Figure 13c. Corresponding, but scaled, areas 532, 534 are blanked off in the 

10 smaller maps. 

Areas larger than the test window may be blanked off in the maps, to avoid 
overlapping detections. In particular, an area equal to the size of the test window surrounded 
by a border half as wide/long as the test window is appropriate to avoid such overlapping 
detections. 

15 Additional faces are detected by searching for the next best response and blanking 

out the corresponding windows successively. 

The intervals allowed between the scales processed are influenced by the sensitivity 
of the method to variations in size. It was found in this preliminary study of scale invariance 
that the method is not excessively sensitive to variations in size as faces which gave a good 

20 response at a certain scale often gave a good response at adjacent scales as well. 

The above description refers to detecting a face even though the size of the face in the 
image is not known at the start of the detection process. Another aspect of multiple scale 
face detection is the use of two or more parallel detections at different scales to validate the 
detection process. This can have advantages if, for example, the face to be detected is 

25 partially obscured, or the person is wearing a hat etc. 

Figures 13d to 13g schematically illustrate this process. During the training phase, 
the system is trained on windows (divided into respective blocks as described above) which 
suiround the whole of the test face (Figure 13d) to generate "full face" histogram data and 
also on windows at an expanded scale so that only a central area of the test face is included 

30 (Figure 13e) to generate "zoomed in" histogram data. This generates two sets of histogram 
data. One set relates to the "full face" windows of Figure 13d, and the other relates to the 
"central face area" windows of Figure 13e. 

During the detection phase, for any given test window 536, the window is applied to 
two different scalings of the test image so that in one (Figure 13f) the test window surrounds 
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the whole of the expected size of a face, and in the other (Figure 13g) the test window 
encompasses the central area of a face at that expected size. These are each processed as 
described above, being compared with the respective sets of histogram data appropriate to 
the type of window. The log probabilities from each parallel process are added before the 
5 comparison with a threshold is applied. 

Putting both of these aspects of multiple scale face detection together leads to a 
particularly elegant saving in the amount of data that needs to be stored. 

In particular, in these embodiments the multiple scales for the arrangements of 
Figures 13a to 13c are arranged in a geometric sequence. In the present example, each scale 

10 in the sequence is a factor of V2 different to the adjacent scale in the sequence. Then, for 
the parallel detection described with reference to Figures 13d to 13g, the larger scale, central 
area, detection is carried out at a scale 3 steps higher in the sequence, that is, 2 % times larger 
than the "full face" scale, using attribute data relating to the scale 3 steps higher in the 
sequence. So, apart from at extremes of the range of multiple scales, the geometric 

15 progression means that the parallel detection of Figures 13d to 13g can always be carried out 
using attribute data generated in respect of another multiple scale three steps higher in the 
sequence. 

The two processes (multiple scale detection and parallel scale detection) can be 
combined in various ways. For example, the multiple scale detection process of Figures 13a 
20 to 13c can be applied first, and then the parallel scale detection process of Figures 13d to 13g 
can be applied at areas (and scales) identified during the multiple scale detection process. 
However, a convenient and efficient use of the attribute data may be achieved by: 

• deriving attributes in respect of the test window at each scale (as in Figures 1 3a to 1 3c) 

• comparing those attributes with the "full face" histogram data to generate a "full face" set 
25 of distance maps 

• comparing the attributes with the "zoomed in" histogram data to generate a "zoomed in" 

set of distance maps 

• for each scale n, combining the "full face" distance map for scale n with the "zoomed in" 
distance map for scale n+3 

30 • deriving face positions from the combined distance maps as described above with 
reference to Figures 13a to 13c 

Further parallel testing can be performed to detect different poses, such as looking 
straight ahead, looking partly up, down, left, right etc. H«re a respective set of histogram 
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data is required and the results are preferably combined using a "max" function, that is, the 
pose giving the highest probability is carried forward to thresholding, the others being 
discarded. 

5 Face Tracking 

A face tracking algorithm will now be described. The tracking algorithm aims to 
improve face detection performance in image sequences. 

The initial aim of the tracking algorithm is to detect every face in every frame of an 
image sequence. However, it is recognised that sometimes a face in the sequence may not be 
10 detected. In these circumstances, the tracking algorithm may assist in interpolating across 
the missing face detections. 

Ultimately, the goal of face tracking is to be able to output some useful metadata 
from each set of frames belonging to the same scene in an image sequence. This might 
include: 
15 • Number of faces. 

• A< Mugshot" (a colloquial word for an image of a person's face, derived from a term 
referring to a police file photograph) of each face. 

• Frame number at which each face first appears. 

• Frame number at which each face last appears. 

20 • Identity of each face (either matched to faces seen in previous scenes, or matched to a 
face database) - this requires some face recognition also. 

The tracking algorithm uses the results of the face detection algorithm, run 
independently on each frame of the image sequence, as its starting point. Because the face 
detection algorithm may sometimes miss (not detect) faces, some method of interpolating the 
25 missing faces is useful. To this end, a K aim an filter was used to predict the next position of 
the face and a skin colour matching algorithm was used to aid tracking of faces. In addition, 
because the face detection algorithm often gives rise to false acceptances, some method of 
rejecting these is also useful. 

The algorithm is shown schematically in Figure 14. 
30 The algorithm will be described in detail below, but in summary, input video data 

545 (representing the image sequence) is supplied to a face detector of the type described in 
this application, and a skin colour matching detector 550. The face detector attempts to 
detect one or more faces in each image. When a face is detected, a Kalman filter 560 is 
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established to track the position of that face. The Kalman filter generates a predicted 
position for the same face in the next image in the sequence. An eye position comparator 
570, 580 detects whether the face detector 540 detects a face at that position (or within a 
certain threshold distance of that position) in the next image. If this is found to be the case, 
then that detected face position is used to update the Kalman filter and the process continues. 

If a face is not detected at or near the predicted position, then a skin colour matching 
method 550 is used. This is a less precise face detection technique which is set up to have a 
lower threshold of acceptance than the face detector 540, so that it is possible for the skin 
colour matching technique to detect (what it considers to be) a face even when the face 
detector cannot make a positive detection at that position. If a "face" is detected by skin 
colour matching, its position is passed to the Kalman filter as an updated position and the 
process continues. 

If no match is found by either the face detector 450 or the skin colour detector 550, 
then the predicted position is used to update the Kalman filter. 

All of these results are subject to acceptance criteria (see below). So, for example, a 
face that is tracked throughout a sequence on the basis of one positive detection and the 
remainder as predictions, or the remainder as skin colour detections, will be rejected. 

A separate Kalman filter is used to track each face in the tracking algorithm. 

In order to use a Kalman filter to track a face, a state model representing the face 
must be created. In the model, the position of each face is represented by a 4-dimensional 
vector containing the co-ordinates of the left and right eyes, which in turn are derived by a 
predetermined relationship to the centre position of the window and the scale being used: 



where k is the frame number. 

The current state of the face is represented by its position, velocity and acceleration, 

in a 12-dimensional vector 



p(k) = 



FirstEyeX 
FirstEyeY 
SecondEyeX 
SecondEyeY 



p(k) 
p(k) 
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First Face Detected 

The tracking algorithm does nothing until it receives a frame with a face detection 
result indicating that there is a face present 

A Kalman filter is then initialised for each detected face in this frame. Its state is 
initialised with the position of the face, and with zero velocity and acceleration: 

>(*)" 



20 



0 
0 



It is also assigned some other attributes: the state model error covariance, Q and the 
observation error covariance, R. The error covariance of the Kalman filter, P, is also 
initialised. These parameters are described in more detail below. At the beginning of the 
10 following frame, and every subsequent frame, a Kalman filter prediction process is carried 
out. 

Kalman Filter Prediction Process 

For each existing Kalman filter, the next position of the face is predicted using the 
15 standard Kalman filter prediction equations shown below. The filter uses the previous state 
(at frame k-1) and some other internal and external variables to estimate the current state of 
the filter (at frame k). 



State prediction equation: z b (k) = <J>(fc, k - \)z a (k - 1) 



Covariance prediction equation: P b (k) = k - l)P a (k - l)o(/r, k - if + Q(k) 

where z b (k) denotes the state before updating the filter for frame k, z a (k-l) denotes the 
state after updating the filter for frame k-l (or the initialised state if it is a new filter), and 
<E>(k 9 k -l) is the state transition matrix. Various state transition matrices were experimented 
25 with, as described below. Similarly, P b (k) denotes the filter's error covariance before 
updating the filter for frame k andP fl (A:-l) denotes the filter's error covariance after 
updating the filter for the previous frame (or the initialised value if it is a new filter). P b (k) 
can be thought of as an internal variable in the filter that models its accuracy. 
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Q{k) is the error covariance of the state model. A high value of Q(k) means that the 
predicted values of the filter's state (i.e. the face's position) will be assumed to have a high 
level of error. By tuning this parameter, the behaviour of the filter can be changed and 
potentially improved for face detection. 



10 



15 



State Transition Matrix 

The state transition matrix, G>(fc,fc-l), determines how the prediction of the next 
state is made. Using the equations for motion, the following matrix can be derived for 
®{k,k-l): 



0(k,k-l) = 



J 4 I A At ±I A {At? 



O a I A 

o A o A 



I A At 



where O a is a 4x4 zero matrix and 7 4 is a 4x4 identity matrix. At can simply be set to 1 (i.e. 

units oft are frame periods). 

This state transition matrix models position, velocity and acceleration. However, it 
was found that the use of acceleration tended to make the face predictions accelerate towards 
the edge of the picture when no face detections were available to correct the predicted state. 
Therefore, a simpler state transition matrix without using acceleration was preferred: 

h 



I A At O a 



o A 

L 0 4 



h 

o A 



o A 
o A 



The predicted eye positions of each Kalman filter, z b {k\ are compared to all face 
20 detection results in the current frame (if there are any). If the distance between the eye 
positions is below a given threshold, then the face detection can be assumed to belong to the 
same face as that being modelled by the Kalman filter. The face detection result is then 
treated as an observation, y(k), of the face's current state: 



><*)' 
0 

0 
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where p{k) is the position of the eyes in the face detection result. This observation is used 
during the Kalman filter update stage to help correct the prediction. 

Skin Colour Matching 

5 Skin colour matching is not used for faces that successfully match face detection 

results. Skin colour matching is only performed for faces whose position has been predicted 
by the Kalman filter but have no matching face detection result in the current frame, and 
therefore no observation data to help update the Kalman filter. 

In a first technique, for each face, an elliptical area centred on the face's previous 

10 position is extracted from the previous frame. An example of such an area 600 within the 
face window 610 is shown schematically in Figure 16. A colour model is seeded using the 
chrominance data from this area to produce an estimate of the mean and covariance of the Cr 
and Cb values, based on a Gaussian model. 

An area around the predicted face position in the current frame is then searched and 

15 the position that best matches the colour model, again averaged over an elliptical area, is 
selected. If the colour match meets a given similarity criterion, then this position is used as 
an observation, y{k) 9 of the face's current state in the same way described for face detection 
results in the previous section. 

Figures 15a and 15b schematically illustrate the generation of the search area. In 

20 particular, Figure 15a schematically illustrates the predicted position 620 of a face within the 
next image 630. In skin colour matching, a search area 640 surrounding the predicted 
position 620 in the next image is searched for the face. 

If the colour match does not meet the similarity criterion, then no reliable observation 
data is available for the current frame. Instead, the predicted state, z b (k) is used as the 

25 observation: 

y{k)=z b {k) 

The skin colour matching methods described above use a simple Gaussian skin 
colour model. The model is seeded on an elliptical area centred on the face in the previous 
frame, and used to find the best matching elliptical area in the current frame. However, to 
30 provide a potentially better performance, two further methods were developed: a colour 
histogram method and a colour mask method. These will now be described. 
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Colour Histogram Method 

In this method, instead of using a Gaussian to model the distribution of colour in the 
tracked face, a colour histogram is used. 

For each tracked face in the previous frame, a histogram of Cr and Cb values within a 
5 square window around the face is computed. To do this, for each pixel the Cr and Cb values 
are first combined into a single value. A histogram is then computed that measures the 
frequency of occurrence of these values in the whole window. Because the number of 
combined Cr and Cb values is large (256x256 possible combinations), the values are 
quantised before the histogram is calculated. 
10 Having calculated a histogram for a tracked face in the previous frame, the histogram 

is used in the current frame to try to estimate the most likely new position of the face by 
finding the area of the image with the most similar colour distribution. As shown 
schematically in Figures 15a and 15b, this is done by calculating a histogram in exactly the 
same way for a range of window positions within a search area of the current frame. This 
15 search area covers a given area around the predicted face position. The histograms are then 
compared by calculating the mean squared error (MSE) between the original histogram for 
the tracked face in the previous frame and each histogram in the current frame. The 
estimated position of the face in the current frame is given by the position of the nfinimum 
MSE. 

20 Various modifications may be made to this algorithm, including: 

• Using three channels (Y, Cr and Cb) instead of two (Cr, Cb). 

• Varying the number of quantisation levels. 

• Dividing the window into blocks and calculating a histogram for each block. In this way, 
the colour histogram method becomes positionally dependent. The MSE between each 

25 pair of histograms is summed in this method. 

• Varying the number of blocks into which the window is divided. 

• Varying the blocks that are actually used - e.g. omitting the outer blocks which might 
only partially contain face pixels. 

For the test data used in empirical trials of these techniques, the best results were 
30 achieved using the following conditions, although other sets of conditions may provide 
equally good or better results with different test data: 

• 3 channels (Y, Cr and Cb). 

• 8 quantisation levels for each channel (i.e. histogram contains 8x8x8 = 512 bins). 
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• Dividing the windows into 16 blocks. 

• Using all 16 blocks. 

Colour Mask Method 

5 This method is based on the method first described above. It uses a Gaussian skin 

colour model to describe the distribution of pixels in the face. 

In the method first described above, an elliptical area centred on the face is used to 
colour match faces, as this may be perceived to reduce or minimise the quantity of 
background pixels which might degrade the model. 

10 In the present colour mask model, a similar elliptical area is still used to seed a colour 

model on the original tracked face in the previous frame, for example by applying the mean 
and covariance of RGB or YCrCb to set parameters of a Gaussian model (or alternatively, a 
default colour model such as a Gaussian model can be used, see below). However, it is not 
used when searching for the best match in the current frame. Instead, a mask area is 

15 calculated based on the distribution of pixels in the original face window from the previous 
frame. The mask is calculated by finding the 50% of pixels in the window which best match 
the colour model. An example is shown in Figures 17a to 17c. In particular, Figure 17a 
schematically illustrates the initial window under test; Figure 17b schematically illustrates 
the elliptical window used to seed the colour model; and Figure 17c schematically illustrates 

20 the mask defined by the 50% of pixels which most closely match the colour model. 

To estimate the position of the face in the current frame, a search area around the 
predicted face position is searched (as before) and the "distance" from the colour model is 
calculated for each pixel. The "distance" refers to a difference from the mean, normalised in 
each dimension by the variance in that dimension. An example of the resultant distance 

25 image is shown in Figure 18. For each position in this distance map (or for a reduced set of 
sampled positions to reduce computation time), the pixels of the distance image are averaged 
over a mask-shaped area. The position with the lowest averaged distance is then selected as 
the best estimate for the position of the face in this frame. 

This method thus differs from the original method in that a mask-shaped area is used 

30 in the distance image, instead of an elliptical area. This allows the colour match method to 
use both colour and shape information. 

Two variations are proposed and were implemented in empirical trials of the 
techniques: 
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(a) Gaussian skin colour model is seeded using the mean and covariance of Cr and Cb 
from an elliptical area centred on the tracked face in the previous frame. 

(b) A default Gaussian skin colour model is used, both to calculate the mask in the 
previous frame and calculate the distance image in the current frame. 
The use of Gaussian skin colour models will now be described further. A Gaussian 

model for the skin colour class is built using the chrominance components of the YCbCr 
colour space. The similarity of test pixels to the skin colour class can then be measured. This 
method thus provides a skin colour likelihood estimate for each pixel, independently of the 
eigenface-based approaches. 
10 Let w be the vector of the CbCr values of a test pixel. The probability of w belonging 

to the skin colour class S is modelled by a two-dimensional Gaussian: 

, , c exp[-^(w-AfcV-^)] 

P(W\S) = —7 

2*|E,|' 

where the mean m s me covariance matrix 2 4 of the distribution are (previously) 

estimated from a training set of skin colour values. 

15 Skin colour detection is not considered to be an effective face detector when used on 

its own. This is because there can be many areas of an image that are similar to skin colour 
but are not necessarily faces, for example other parts of the body. However, it can be used to 
improve the performance of the eigenblock-based approaches by using a combined approach 
as described in respect of the present face tracking system. The decisions made on whether 

20 to accept the face detected eye positions or the colour matched eye positions as the 
observation for the Kalman filter, or whether no observation was accepted, are stored. These 
are used later to assess the ongoing validity of the faces modelled by each Kalman filter. 

Kalman Filter Update Step 
25 The update step is used to determine an appropriate output of the filter for the current 

frame, based on the state prediction and the observation data. It also updates the internal 
variables of the filter based on the error between the predicted state and the observed state. 
The following equations are used in the update step: 
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Kalman gain equation 
State update equation 



K(k) = P„ {k)H T {kiH{k)P b (k)H T (*) + R{k)Y 
PM = P b (k)-K(k)H(k)P b (k) 



Covariance update equation 

Here, K(k) denotes the Kalman gain, another variable internal to the Kalman filter. It 
is used to determine how much the predicted state should be adjusted based on the observed 
state, y{k). 

H(k) is the observation matrix. It determines which parts of the state can be 
observed. In our case, only the position of the face can be observed, not its velocity or 
acceleration, so the following matrix is used for H (k) : 



H(k) = 



h o A o 4 
o 4 o 4 o 4 
o 4 o 4 o 4 



25 



R(k) is the error covariance of the observation data. In a similar way to Q(k), a high 
value of R(k) means that the observed values of the filter's state (i.e. the face detection 
results or colour matches) will be assumed to have a high level of error. By tuning this 
parameter, the behaviour of the filter can be changed and potentially improved for face 
detection. For our experiments, a large value of R(k) relative to Q(k) was found to be 
suitable (this means that the predicted face positions are treated as more reliable than the 
observations). Note that it is permissible to vary these parameters from frame to frame. 
Therefore, an interesting future area of investigation may be to adjust the relative values of 
R(k) and Q(k) depending on whether the observation is based on a face detection result 
(reliable) or a colour match (less reliable). 

For each Kalman filter, the updated state, z a {k), is used as the final decision on the 

position of the face. This data is output to file and stored. 

Unmatched face detection results are treated as new faces. A new Kalman filter is 
initialised for each of these. Faces are removed which: 
• Leave the edge of the picture and/or 
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. Have a lack of ongoing evidence supporting them (when there is a high proportion of 
observations based on Kalman filter predictions rather than face detection results or 
colour matches). 

For these faces, the associated Kalman filter is removed and no data is output to file. 
As an optional difference from this approach, where a face is detected to leave the picture, 
the tracking results up to the frame before it leaves the picture may be stored and treated as 
valid face tracking results (providing that the results meet any other criteria applied to 

validate tracking results). 

These rules may be formalised and built upon by bringing in some additional 

variables: 



predictionjxcceptancejratiojhreshold 



If, during tracking a given face, the proportion 
of accepted Kalman predicted face positions 
exceeds this threshold, then the tracked face is 
rejected. 

This is currently set to 0.8. 



20 



detection jicceptancejraUo ^threshold 



During a final pass through all the frames, if for 
a given face the proportion of accepted face 
detections falls below this threshold, then the 
tracked face is rejected. 
This is currently set to 0.08. 



min_ frames 



25 



During a final pass through all the frames, if for 
a given face the number of occurrences is less 
than min_frames, the face is rejected. This is 
only likely to occur near the end of a sequence, 
minframes is currently set to 5. 



30 final^rediction^acceptance^atiojhreshold and minjrames2 During a final pass 

through all the frames, if for a given tracked 
face the number of occurrences is less than 
min_frames2 AND the proportion of accepted 
Kalman predicted face positions exceeds the 
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final jrediction_acceptance_ratio_threshold, 
the face is rejected. Again, this is only likely to 
occur near the end of a sequence. 
finaljprediction_accept^ce_ratio_threshold is 
currently set to 0.5 and min_frames2 is 
currently set to 10. 

minjeye_spacing Additionally, faces are now removed if they are 

tracked such that the eye spacing is decreased 
below a given minimum distance. This can 
happen if the Kalman filter falsely believes the 
eye distance is becoming smaller and there is no 
other evidence, e.g. face detection results, to 
correct this assumption. If uncorrected, the eye 
distance would eventually become zero. As an 
optional alternative, a minimum or lower limit 
eye separation can be forced, so that if the 
detected eye separation reduces to the minimum 
eye separation, the detection process continues 
to search for faces having that eye separation, 
but not a smaller eye separation. 

It is noted that the tracking process is not limited to tracking through a video 
sequence in a forward temporal direction. Assuming that the image data remain accessible 
(i.e. the process is not real-time, or the image data are buffered for temporary continued use), 
the entire tracking process could be carried out in a reverse temporal direction. Or, when a 
first face detection is made (often part-way through a video sequence) the tracking process 
could be initiated in both temporal directions. As a further option, the tracking process could 
be run in both temporal directions through a video sequence, with the results being combined 
so that (for example) a tracked face meeting the acceptance criteria is included as a valid 
result whichever direction the tracking took place. 
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Overlap Rules for Face Tracking 

When the faces are tracked, it is possible for the face tracks to become overlapped. 
When this happens, in at least some applications, one of the tracks should be deleted. A set 
of rules is used to determine which face track should persist in the event of an overlap. 

Whilst the faces are being tracked there are 3 possible types of track: 

D: Face Detection - the current position of the face is confirmed by a new face 
detection 

S: Skin colour track - there is no face detection, but a suitable skin colour track 
has been found 

P: Prediction - there is neither a suitable face detection nor skin colour track, so 
the predicted face position from the Kalman filter is used. 

The following grid defines a priority order if two face tracks overlap with each other: 



~~ Face 2 
Face 1 ' -- 


D 


S 


P 


D 


Largest Face Size 


D 


D 


S 


D 


Largest Face Size 


S 


P 


D 


S 


Largest Face Size 



So, if both tracks are of the same type, then the largest face size determines which 
track is to be maintained. Otherwise, detected tracks have priority over skin colour or 
predicted tracks. Skin colour tracks have priority over predicted tracks. 

15 In the tracking method described above, a face track is started for every face 

detection that cannot be matched up with an existing track. This could lead to many false 
detections being erroneously tracked and persisting for several frames before finally being 
rejected by one of the existing rules (e.g. by the rule associated with the 
prediction_aeceptwice_ratio_threskold) 

20 Also, the existing rules for rejecting a track (e.g. those rules relating to the variables 

prediction_acceptance_ratioJhreshold and ,detection_acceptance_ratio_threshold), are 
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biased against tracking someone who turns their head to the side for a significant length of 
time. In reality, it is often desirable to carry on tracking someone who does this. 
A solution will now be described. 

The first part of the solution helps to prevent false detections from setting off 
erroneous tracks. A face track is still started internally for every face detection that does not 
match an existing track. However, it is not output from the algorithm. In order for this track 
to be maintained, the first / frames in the track must be face detections (i.e. of type D ). If all 
of the first /frames are of type D then the track is maintained and face locations are output 
from the algorithm from frame/ onwards. 

If all of the first n frames are not of type D, then the face track is terminated and no 
face locations are output for this track. 

/ is typically set to 2, 3 or 5. 

The second part of the solution allows faces in profile to be tracked for a long period, 
rather than having their tracks terminated due to a low detection_acceptance_ratio. To 
achieve this, where the faces are matched by the ± 30° eigenblocks, the tests relating to the 
variables prediction jicceptance^ratiojthreshold and detection_acceptance_ratiojthreshdld 
are not used. Instead, an option is to include the following criterion to maintain a face track: 
g consecutive face detections are required every n frames to maintain the face track 
where g is typically set to a similar value tof 9 e.g. 1-5 frames and n corresponds to 
the maximum number of frames for which we wish to be able to track someone when they 
are turned away from the camera, e.g. 10 seconds (= 250 or 300 frames depending on frame 
rate). 

This may also be combined with the prediction _acceptance^ratio_threshold and 
detection_acceptance_ratio_jhreshold rules. Alternatively, the 

predictionjxcceptancejratiojhreshold and detection jxcceptancej-atiojhreshold may be 
applied on a rolling basis e.g. over only the last 30 frames, rather than since the beginning of 
the track. 

Another criterion for rejecting a face track is that a so-called "bad colour threshold" 
is exceeded. In this test a tracked face position is validated by skin colour (whatever the 
acceptance type - face detection or Kalman prediction). Any face whose distance from an 
expected skin colour exceeds a given "bad colour threshold" has its track terminated. 

In the method described above, the skin colour of the face is only checked during 
skin colour tracking. This means that non-skin-coloured false detections may be tracked, or 
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the face track may wander off into non-skin-coloured locations by using the predicted face 
position. 

To improve on mis, whatever the acceptance type of the face (detection, skin colour 
or Kalman prediction), its skin colour is checked. If its distance (difference) from skin colour 
5 exceeds a bad_colour_threshold, then the face track is terminated. 

An efficient way to implement this is to use the distance from skin colour of each 
pixel calculated during skin colour tracking. If this measure, averaged over the face area 
(either over a mask shaped area, over an elliptical area or over the whole face window 
depending on which skin colour tracking method is being used), exceeds a fixed threshold, 
1 o then the face track is terminated. 

A further criterion for rejecting a face track is that its variance is very low or very 
high. This technique will be described below after the description of Figures 22a to 22c. 

In the tracking system shown schematically in Figure 14, three further features are 
included. 

Shot boundary data 560 (from metadata associated with the image sequence under 
test; or metadata generated within the camera of Figure 2) defines the limits of each 
contiguous "shot" within the image sequence. The Kalman filter is reset at shot boundaries, 
and is not allowed to carry a prediction over to a subsequent shot, as the prediction would be 
meaningless. 

User metadata 542 and camera setting metadata 544 are supplied as inputs to the face 
detector 540: These may also be used in a non-tracking system. Examples of the camera 
setting metadata were described above. User metadata may include information such as: 

• type of programme (e.g. news, interview, drama) 

• script information such as specification of a "long shot" , "medium close-up" etc 
(particular types of camera shot leading to an expected sub-range of face sizes), how 
many people involved in each shot (again leading to an expected sub-range of face sizes) 
and so on 

• sports-related information - sports are often filmed from fixed camera positions using 
standard views and shots. By specifying these in the metadata, again a sub-range of face 

30 sizes can be derived 

The type of programme is relevant to the type of face which may be expected in the 
images or image sequence. For example, in a news programme, one would expect to see a 
single face for much of the image sequence, occupying an area of (say) 10% of the screen. 
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The detection of faces at different scales can be weighted in response to this data, so that 
faces of about this size are given an enhanced probability. Another alternative or additional 
approach is that the search range is reduced, so that instead of searching for faces at all 
possible scales, only a subset of scales is searched. This can reduce the processing 
requirements of the face detection process. In a software-based system, the software can run 
more quickly and/or on a less powerful processor. In a hardware-based system (including 
for example an application-specific integrated circuit (ASIC) or field programmable gate 
array (FPGA) system) the hardware needs may be reduced. 

The other types of user metadata mentioned above may also be applied in this way. 
The "expected face size" sub-ranges may be stored in a look-up table held in the memory 30, 
for example. 

As regards camera metadata, for example the current focus and zoom settings of the 
lens 110, these can also assist the face detector by giving an initial indication of the expected 
image size of any faces that may be present in the foreground of the image. In this regard, it 
is noted that the focus and zoom settings between them define the expected separation 
between the camcorder 100 and a person being filmed, and also the magnification of the lens 
110. From these two attributes, based upon an average face size, it is possible to calculate 
the expected size (in pixels) of a face in the resulting image data, leading again to a sub- 
range of sizes for search or a weighting of the expected face sizes. 

This arrangement lends itself to use in a video conferencing or so-called digital 
signage environment. 

In a video conferencing arrangement the user could classify the video material as 
"individual speaker", "Group of two", "Group of three" etc, and based on this classification a 
face detector could derive an expected face size and could search for and highlight the one or 
more faces in the image. 

In a digital signage environment, advertising material could be displayed on a video 
screen. Face detection is used to detect the faces of people looking at the advertising 
material. 

Advantages of the tracking algorithm 

The face tracking technique has three main benefits: 
• It allows missed faces to be filled in by using Kalman filtering and skin colour tracking 

in frames for which no face detection results are available. This increases the true 

acceptance rate across the image sequence. 
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• It provides face linking: by successfully tracking a face, the algorithm automatically 
knows whether a face detected in a future frame belongs to the same person or a different 
person. Thus, scene metadata can easily be generated from this algorithm, comprising the 
number of faces in the scene, the frames for which they are present and providing a 

5 representative mugshot of each face. 

• False face detections tend to be rejected, as such detections tend not to carry forward 
between images. 

Figures 19a to 19c schematically illustrate the use of face tracking when applied to a 
video scene. 

10 In particular, Figure 19a schematically illustrates a video scene 800 comprising 

successive video images (e.g. fields or frames) 810. 

In this example, the images 810 contain one or more faces. In particular all of the 

images 810 in the scene include a face A, shown at an upper left-hand position within the 

schematic representation of the image 810. Also, some of the images include a face B shown 
15 schematically at a lower right hand position within the schematic representations of the 

images 810. 

A face tracking process is applied to the scene of Figure 19a. Face A is tracked 
reasonably successfully throughout the scene. In one image 820 the face is not tracked by a 
direct detection, but the skin colour matching techniques and the Kalman filtering techniques 

20 described above mean that the detection can be continuous either side of the "missing" image 
820. The representation of Figure 19b indicates the detected probability of a face being 
present in each of the images. It can be seen that the probability is highest at an image 830, 
and so the part 840 of the image detected to contain face A is used as a "picture stamp" in 
respect of face A. Picture stamps will be described in more detail below. 

25 Similarly, face B is detected with different levels of confidence, but an image 850 

gives rise to the highest detected probability of face B being present. Accordingly, the part 
of the corresponding image detected to contain face B (part 860) is used as a picture stamp 
for face B within that scene. (Alternatively, of course, a wider section of the image, or even 
the whole image, could be used as the picture stamp). 

30 For each tracked face, a single representative face picture stamp is required. 

Outputting the face picture stamp based purely on face probability does not always give the 
best quality of picture stamp. To get the best picture quality it would be better to bias or 
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steer the selection decision towards faces that are detected at the same resolution as the 

picture stamp, e.g. 64x64 pixels 

To get the best quality picture stamps the following scheme may be applied: 
(1) Use a face that was detected (not colour tracked / Kalman tracked) 
5 (2) Use a face that gave a high probability during face detection, i.e. at least a 

threshold probability 

(3) Use a face which is as close as possible to 64x64 pixels, to reduce rescaling 
artefacts and improve picture quality 

(4) Do not (if possible) use a very early face in the track, i.e. a face in a 
10 predetermined initial portion of the tracked sequence (e.g. 10% of the tracked sequence, or 

20 frames, etc) in case this means that the face is still very distant (i.e. small) and blurry 
Some rules that could achieve this are as follows: 
For each face detection: 

15 Calculate the metric M = face_probability * size_weighting, where size_weighting = 

MIN( (face_size/64y% (64/face_size) A x) and x=0.25. Then take the face picture stamp for 
which M is largest. 

This gives the following weightings on the face probability for each face size: 

20 



facesize 


size_weighting 


16 


0.71 


19 


0.74 


23 


0.77 


27 


0.81 


32 


0.84 


38 


0.88 


45 


0.92 


54 


0.96 


64 


1.00 


76 


0.96 


91 


0.92 


108 


0.88 


128 


0.84 
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152 


0.81 


181 


0.77 


215 


0.74 


256 


0.71 


304 


0.68 


362 


0.65 


431 


0.62 


512 


0.59 



10 In practice this could be done using a look-up table. 

To make the weighting function less harsh, a smaller power than 0.25, e.g x=0.2 or 
0.1, could be used. 

This weighting technique could be applied to the whole face track or just to the first 
15 N frames (to apply a weighting against the selection of a poorly-sized face from those N 
frames). N could for example represent just the first one or two seconds (25-50 frames). 

In addition, preference is given to faces that are frontally detected over those that 
were detected at +- 30 degrees (or any other pose). 

Figure 20 schematically illustrates a display screen of a non-linear editing system. 
20 Non-linear editing systems are well established and are generally implemented as 

software programs running on general purpose computing systems such as the system of 
Figure 1. These editing systems allow video, audio , and other material to be edited to an 
output media product in a manner which does not depend on the order in which the 
individual media items (e.g. video shots) were captured. 
25 The schematic display screen of Figure 20 includes a viewer area 900, in which video 

clips be may viewed, a set of clip icons 910, to be described further below and a * timeline" 
920 including representations of edited video shots 930, each shot optionally containing a 
picture stamp 940 indicative of the content of that shot. 

At one level, the face picture stamps derived as described with reference to Figures 
30 19a to 19c could be used as the picture stamps 940 of each edited shot so, within the edited 
length of the shot, which may be shorter than the originally captured shot, the picture stamp 
representing a face detection which resulted in the highest face probability value can be 
inserted onto the time line to show a representative image from that shot. The probability 
values may be compared with a threshold, possibly higher than 4he basic face detection 
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threshold, so that only face detections having a high level of confidence are used to generate 
picture stamps in this way. If more than one face is detected in the edited shot, the face with 
the highest probability may be displayed, or alternatively more than one face picture stamp 
may be displayed on the time line. 

5 Time lines in non-linear editing systems are usually capable of being scaled, so that 

the length of line corresponding to the full width of the display screen can represent various 
different time periods in the output media product So, for example, if a particular boundary 
between two adjacent shots is being edited to frame accuracy, the time line may be 
"expanded" so that the width of the display screen represents a relatively short time period in 

10 the output media product. On the other hand, for other purposes such as visualising an 
overview of the output media product, the time line scale may be contracted so that a longer 
time period may be viewed across the width of the display screen. So, depending on the 
level of expansion or contraction of the time line scale, there may be less or more screen area 
available to display each edited shot contributing to the output media product. 

15 In an expanded time line scale, there may well be more than enough room to fit one 

picture stamp (derived as shown in Figures 19a to 19c) for each edited shot making up the 
output media product. However, as the time line scale is contracted, this may no longer be 
possible. In such cases, the shots may be grouped together in to "sequences", where each 
sequence is such that it is displayed at a display screen size large enough to accommodate a 

20 phase picture stamp. From within the sequence, then, the face picture stamp having the 
highest corresponding probability value is selected for display. If no face is detected within 
a sequence, an arbitrary image, or no image, can be displayed on the timeline. 

Figure 20 also shows schematically two "face timelines" 925, 935. These scale with 
the "main" timeline 920. Each face timeline relates to a single tracked face, and shows the 

25 portions of the output edited sequence containing that tracked face. It is possible that the 
user may observe that certain faces relate to the same person but have not been associated 
with one another by the tracking algorithm. The user can "link" these faces by selecting the 
relevant parts of the face timelines (using a standard Windows*™ selection technique for 
multiple items) and then clicking on a "link" screen button (not shown). The face timelines 

30 would then reflect the linkage of the whole group of face detections into one longer tracked 
face. Figures 21a and 21b schematically illustrate two variants of clip icons 910' and 910". 
These are displayed on the display screen of Figure 20 to allow the user to select individual 
clips for inclusion in the time line and editing of their start and end positions (in and out 
points). So, each clip icon represents the whole of a respective clip stored on the system. 
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hi Figure 21a, a clip icon 910" is represented by a single face picture stamp 912 and 
a text label area 914 which may include, for example, time code information defining the 
position and length of that clip. In an alternative arrangement shown in Figure 21b, more 
than one face picture stamp 916 may be included by using a multi-part clip icon. 
5 Another possibility for the clip icons 910 is that they provide a "face summary" so 

that all detected faces are shown as a set of clip icons 910, in the order in which they appear 
(either in the source material or in the edited output sequence). Again, faces that are the 
same person but which have not been associated with one another by the tracking algorithm 
can be linked by the user subjectively observing that they are the same face. The user could 
10 select the relevant face clip icons 910 (using a standard Windows*™ selection technique for 
multiple items) and then click on a "link" screen button {not shown). The tracking data 
would then reflect the linkage of the whole group of face detections into one longer tracked 
face. 

A further possibility is that the clip icons 910 could provide a hyperlink so that the 
15 user may click on one of the icons 910 which would then cause the corresponding clip to be 
played in the viewer area 900. 

A similar technique may be used in, for example, a surveillance or closed circuit 
television (CCTV) system. Whenever a face is tracked, or whenever a face is tracked for at 
least a predetermined number of frames, an icon similar to a clip icon 910 is generated in 
20 respect of the continuous portion of video over which that face was tracked. The icon is 
displayed in a similar manner to the clip icons in Figure 20. Clicking on an icon causes the 
replay (in a window similar to the viewer area 900) of the portion of video over which that 
particular face was tracked. It will be appreciated that multiple different faces could be 
tracked in this way, and that the corresponding portions of video could overlap or even 

25 completely coincide. 

Figures 22a to 22c schematically illustrate a gradient pre-processing technique. 

It has been noted that image windows showing little pixel variation ^can tend to be 
detected as faces by a face detection arrangement based on eigenfaces or eigeriblocks. 
Therefore, a pre-processing step is proposed to remove areas of little pixel variation from the 
30 face detection process. In the case of a multiple scale system (see above) the pre-processing 
step can be carried out at each scale. 

The basic process is that a "gradient test" is applied to each possible window position 
across the whole image. A predetermined pixel position for each window position, such as 
the pixel at or nearest the centre of that window position, is flagged or labelled in 
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dependence on the results of the test applied to that window. If the test shows that a window 
has little pixel variation, that window position is not used in the face detection process. 

A first step is illustrated in Figure 22a. This shows a window at an arbitrary window 
position in the image. As mentioned above, the pre-processing is repeated at each possible 

5 window position. Referring to Figure 22a, although the gradient pre-processing could be 
applied to the whole window, it has been found that better results are obtained if the pre- 
processing is applied to a central area 1000 of the test window 1010. 

Referring to Figure 22b, a gradient-based measure is derived from the window (or 
from the central area of the window as shown in Figure 22a), which is the average of the 

10 absolute differences between all adjacent pixels 1011 in both the horizontal and vertical 
directions, taken over the window. Each window centre position is labelled with this 
gradient-based measure to produce a gradient "map" of the image. The resulting gradient 
map is then compared with a threshold gradient value. Any window positions for which the 
gradient-based measure lies below the threshold gradient value are excluded from the face 

15 detection process in respect of that image. 

Alternative gradient-based measures could be used, such as the pixel variance or the 
mean absolute pixel difference from a mean pixel value. 

The gradient-based measure is preferably carried out in respect of pixel luminance 
values, but could of course be applied to other image components of a colour image. 

20 Figure 22c schematically illustrates a gradient map derived from an example image. 

Here a lower gradient area 1070 (shown shaded) is excluded from face detection, and only a 
higher gradient area 1080 is used. The embodiments described above have related to a 
face detection system (involving training and detection phases) and possible uses for it in a 
camera-recorder and an editing system. It will be appreciated that there are many other 

25 possible uses of such techniques, for example (and not limited to) security surveillance 
systems, media handling in general (such as video tape recorder controllers), video 
conferencing systems and the like. 

In other embodiments, window positions having high pixel differences can also be 
flagged or labelled, and are also excluded from the face detection process. A **high" pixel 

30 difference means that the measure described above with respect to Figure 22b exceeds an 
upper threshold value. 

So, a gradient map is produced as described above. Any positions for which the 
gradient measure is lower than the (first) threshold gradient value mentioned earlier are 
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excluded from face detection processing, as are any positions for which the gradient measure 
is higher than the upper threshold value. 

It was mentioned above that the "lower threshold" processing is preferably applied to 
a central part 1000 of the test window 1010. The same can apply to the "upper threshold" 

5 processing. This would mean that only a single gradient measure needs to be derived in 
respect of each window position. Alternatively, if the whole window is used in respect of the 
lower threshold test, the whole window can similarly be used in respect of the upper 
threshold test. Again, only a single gradient measure needs to be derived for each window 
position. Of course, however, it is possible to use two different arrangements, so that (for 

10 example) a central part 1000 of the test window 1010 is used to derive the gradient measure 
for the lower threshold test, but the full test window is used in respect of the upper threshold 
test. 

A further criterion for rejecting a face track, mentioned earlier, is mat its variance or 
gradient measure is very low or very high. 
15 In this technique a tracked face position is validated by variance from area of interest 

map. Only a face-sized area of the map at the detected scale is stored per face for the next 
iteration of tracking. 

Despite the gradient pre-processing described above, it is still possible for a skin 
colour tracked or Kalman predicted face to move into a (non-face-like) low or high variance 
20 area of the image. So, during gradient pre-processing, the variance values (or gradient 
values) for the areas around existing face tracks are stored. 

When the final decision on the face's next position is made (with any acceptance 
type, either face detection, skin colour or Kalman prediction) the position is validated against 
the stored variance (or gradient) values in the area of interest map. If the position is found to 
25 have very high or very low variance (or gradient), it is considered to be non-face-like and the 
face track is terminated. This prevents face tracks from wandering onto low (or high) 
variance background areas of the image. 

Alternatively, even if gradient pre-processing is not used, the variance of the new 
face position can be calculated afresh. In either case the variance measure used can either be 
30 traditional variance or the sum of differences of neighbouring pixels (gradient) or any other 
variance-type measure. 

Figure 23 schematically illustrates a video conferencing system. Two video 
conferencing stations 1100, 1110 are connected by a network connection 1120 such as: the 
Internet, a local or wide area network, a telephone line, a high bit rate leased line, an ISDN 
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line etc. Each of the stations comprises, in simple terms, a camera and associated sending 
apparatus 1130 and a display and associated receiving apparatus 1140. Participants in the 
video conference are viewed by the camera at their respective station and their voices are 
picked up by one or more microphones (not shown in Figure 23) at that station. The audio 
and video information is transmitted via the network 1120 to the receiver 1140 at the other 
station. Here, images captured by the camera are displayed and the participants' voices are 
produced on a loudspeaker or the like. 

It will be appreciated that more than two stations may be involved in the video 
conference, although the discussion here will be limited to two stations for simplicity. 

Figure 24 schematically illustrates one channel, being the connection of one 
camera/sending apparatus to one display/receiving apparatus. 

At the camera/sending apparatus, there is provided a video camera 1150, a face 
detector 1160 using the techniques described above, an image processor 1170 and a data 
formatter and transmitter 1180. A microphone 1 190 detects the participants' voices. 

Audio, video and (optionally) metadata signals are transmitted from the formatter and 
transmitter 1180, via the network connection 1120 to the display/receiving apparatus 1140. 
Optionally, control signals are received via the network connection 1120 from the 
display/receiving apparatus 1 140. 

At the display/receiving apparatus, there is provided a display and display processor 
1200, for example a display screen and associated electronics, user controls 1210 and an 
audio output arrangement 1220 such as a digital to analogue (D AC) converter, an amplifier 
and a loudspeaker. 

In general terms, the face detector 1160 detects (and optionally tracks) faces in the 
captured images from the camera 1150. The face detections are passed as control signals to 
the image processor 1 170. The image processor can act in various different ways, which will 
be described below, but fundamentally the image processor 1170 alters the images captured 
by the camera 1 150 before they are transmitted via the network 1 120. A significant purpose 
behind this is to make better use of the available bandwidth or bit rate which can be carried 
by the network connection 1 120. Here it is noted that in most commercial applications, the 
cost of a network connection 1 120 suitable for video conference purposes increases with an 
increasing bit rate requirement. At the formatter and transmitter 1180 the images from the 
image processor 1170 are combined with audio signals from the microphone 1190 (for 
example, having been converted via an analogue to digital converter (ADC)) and optionally 
metadata defining the nature of the processing carried out by the image processor 1 170. 
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Various modes of operation of the video conferencing system will be described 

below. 

Figure 25 is a further schematic representation of the video conferencing system. 
Here, the functionality of the face detector 1 160, the image processor 1 170, the formatter and 
5 transmitter 1180 and the processor aspects of the display and display processor 1200 are 
carried out by programmable personal computers 1230. The schematic displays shown on 
the display screens (part of 1200) represent one possible mode of video conferencing using 
face detection which will be described below with reference to Figure 31, namely that only 
those image portions containing faces are transmitted from one location to the other, and are 
10 then displayed in a tiled or mosaic form at the other location. As mentioned, this mode of 
operation will be discussed below. 

Figure 26 is a flowchart schematically illustrating a mode of operation of the system 
of Figures 23 to 25. The flowcharts of Figures 26, 28, 31, 33 and 34 are divided into 
operations carried out at the camera/sender end (1130) and those carried out at the 
15 display/receiver end (1 140). 

So, referring to Figure 26, the camera 1 150 captures images at a step 1300. At a step 
1310, the face detector 1 160 detects faces in the captured images. Ideally, face tracking (as 
described above) is used to avoid any spurious interruptions in the face detection and to 
provide that a particular person's face is treated in the same way throughout the video 
20 conferencing session. 

At a step 1320, the image processor 1170 crops the captured images in response to 
the face detection information. This may be done as follows: 

first, identify the upper left-most face detected by the face detector 1 160 
detect the upper left-most extreme of that face; this forms the upper left 
25 comer of the cropped image 

repeat for the lower right-most face and the lower right-most extreme of that 
face to form the lower right comer of the cropped image 

crop the image in a rectangular shape based on these two co-ordinates. 

30 The cropped image is then transmitted by the formatter and transmitted 1 1 80. In this 

instance, there is no need to transmit additional metadata. The cropping of the image allows 
either a reduction in bit rate compared to the full image or an improvement in transmission 
quality while maintaining the same bit rate. 
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At the receiver, the cropped image is displayed at a full screen display at a step 1 130. 

Optionally, a user control 1210 can toggle the image processor 1 170 between a mode 
in which the image is cropped and a mode in which it is not cropped. This can allow the 
participants at the receiver end to see either the whole room or just the face-related parts of 
5 the image. 

Another technique for cropping the image is as follows: 

• identify the leftmost and rightmost faces 

• maintaining the aspect ratio of the shot, locate the faces in the upper half of 
the picture. 

10 In an alternative to cropping, the camera could be zoomed so that the detected faces 

are featured more significantly in the transmitted images. This could, for example, be 
combined with a bit rate reduction technique on the resulting image. To achieve this, a 
control of the directional (pan/tilt) and lens zoom properties of the camera, is made available 
to the image processor (represented by a dotted line 1 155 in Figure 24) 

15 Figures 27a and 27b are example images relating to the flowchart of Figure 26. 

Figure 27a represents a full screen image as captured by the camera 1150, whereas Figure 
27b represents a zoomed version of that image. 

Figure 28 is a flowchart schematically illustrating another mode of operation of the 
system of Figures 23 to 25. Step 1300 is the same as that shown in Figure 26. 

20 At a step 1340, each face in the captured images is identified and highlighted, for 

example by drawing a box around that face for display. Each face is also labelled, for 
example with an arbitrary label a, b, c... Here, face tracking is particularly useful to avoid 
any subsequent confusion over the labels. The labelled image is formatted and transmitted to 
the receiver where it is displayed at a step 1350. At a step 1360, the user selects a face to be 

25 displayed, for example by typing the label relating to that face. The selection is passed as 
control data back to the image processor 1 170 which isolates the required face at a step 1370. 
The required face is transmitted to the receiver. At a step 1380 the required face is 
displayed. The user is able to select a different face by the step 1360 to replace the currently 
displayed face. Again, this arrangement allows a potential saving in bandwidth, in that the 

30 selection screen may be transmitted at a lower bit rate because it is only used for selecting a 
face to be displayed. Alternatively, as before, the individual faces, once selected, can be 
transmitted at an enhanced bit rate to achieve a better quality image. 
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Figure 29 is an example image relating to the flowchart of Figure 28. Here, three 
faces have been identified, and are labelled a, b and c. By typing one of those three letters 
into the user controls 1210, the user can select one of those faces for a full-screen display. 
This can be achieved by a cropping of the main image or by the camera zooming onto that 
5 face as described above. Figure 30 shows an alternative representation, in which so-called 
thumbnail images of each face are displayed as a menu for selection at the receiver. 

Figure 31 is a flowchart schematically illustrating a further mode of operation of the 
system of Figures 23 to 25. The steps 1300 and 1310 correspond to those of Figure 26. 

At a step 1400, the image processor 1170 and the formatter and transmitter 1180 co- 
10 operate to transmit only thumbnail images relating to the captured faces. These are 
displayed as a menu or mosaic of faces at the receiver end at a step 1410. At a step 1420, 
optionally, the user can select just one face for enlarged display. This may involve keeping 
the other faces displayed in a smaller format on the same screen or the other faces may be 
hidden while the enlarged display is used. So a difference between this arrangement and that 
15 of Figure 28 is that thumbnail images relating to all of the faces are transmitted to the 
receiver, and the selection is made at the receiver end as to how the thumbnails are to be 
displayed. 

Figure 32 is an example image relating to the flowchart of Figure 31. Here, an initial 
screen could show three thumbnails, 1430, but the stage illustrated by Figure 32 is that the 
20 face belonging to participant c has been selected for enlarged display on a left hand part of 
the display screen. However, the thumbnails relating to the other participants are retained so 
that the user can make a sensible selection of a next face to be displayed in enlarged form. 

It should be noted that, at least in a system where the main image is cropped, the 
thumbnail images referred to in these examples are "live" thumbnail images, albeit taking 
25 into account any processing delays present in the system. That is to say, the thumbnail 
images vary in time, as the captured images of the participants vary. In a system using a 
camera zoom, then the thumbnails could be static or a second camera could be used to 
capture the wider angle scene. 

Figure 33 is a flowchart schematically illustrating a further mode of operation. Here, 
30 the steps 1300 and 1310 correspond to those of Figure 26. 

At a step 1440 a thumbnail face image relating to the face detected to be nearest to an 
active microphone is transmitted. Of course, this relies on having more than one microphone 
and also a pre-selection or metadata defining which participant is sitting near to which 
microphone. This can be set up in advance by a simple menu-driven table entry by the users 
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at each video conferencing station. The active microphone is considered to be the 
microphone having the greatest magnitude audio signal averaged over a certain time (such as 
one second). A low pass filtering arrangement can be used to avoid changing the active 
microphone too often, for example in response to a cough or an object being dropped, or two 
5 participants speaking at the same time. 

At a step 1450 the transmitted face is displayed. A step 1460 represents the quasi- 
continuous detection of a current active microphone. 

The detection could be, for example, a detection of a single active microphone or 
alternatively a simple triangulation technique could detect the speaker's position based on 
10 multiple microphones. 

Finally, Figure 34 is a flowchart schematically illustrating another mode of operation, 
again in which the steps 1300 and 1310 correspond to those of Figure 26. 

At a step 1470 the parts of the captured images immediately surrounding each face 
are transmitted at a higher resolution and the background (other parts of the captured images) 
15 is transmitted at a lower resolution. This can achieve a useful saving in bit rate or allow an 
enhancement of the parts of the image surrounding each face. Optionally, metadata can be 
transmitted defining the position of each face, or the positions may be derived at the receiver 
by noting the resolution of different parts of the image. 

At a step 1480, at the receiver end the image is displayed and the faces are optionally 
20 labelled for selection by a user at a step 1490 this selection could cause the selected face to 
be displayed in a larger format similar to the arrangement of Figure 32. 

Although the description of Figures 23 to 34 has related to video conferencing 
systems, the same techniques could be applied to, for example, security monitoring (CCTV) 
systems. Here, a return channel is not normally required, but an arrangement as shown in 
25 Figure 24, where the camera / sender arrangement is provided as a CCTV camera, and the 
receiver / display arrangement is provided at a monitoring site, could use the same 
techniques as those described for video conferencing. 

It will be appreciated that the embodiments of the invention described above may of 
course be implemented, at least in part, using software-controlled data processing apparatus. 
30 For example, one or more of the components schematically illustrated or described above may 
be implemented as a software-controlled general purpose data processing device or a bespoke 
program controlled data processing device such as an application specific integrated circuit, a 
field programmable gate array or the like. It will be appreciated that a computer program 
providing such software or program control and a storage, transmission or other providing 
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medium by which such a computer program is stored are envisaged as aspects of the present 
invention. 

The list of references and appendices follow. For the avoidance of doubt, it is noted 
5 that the list and the appendices form a part of the present description. These documents are 
all incorporated by reference. 
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Appendix A: Training Face Sets 

One database consists of many thousand images of subjects standing in front of an indoor 
25 background. Another training database used in experimental implementations of the above 
techniques consists of more than ten thousand eight-bit greyscale images of human heads 
with views ranging from frontal to left and right profiles. The skilled man will of course 
understand that various different training sets could be used, optionally being profiled to 
reflect facial characteristics of a local population. 
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Appendix B - Eigenblocks 

In the eigenface approach to face detection and recognition (References 4 and 5), 

5 each /w-by-n face image is reordered so that it is represented by a vector of length mn. Each 
image can then be thought of as a point in wm-dimensional space. A set of images maps to a 
collection of points in this large space. 

Face images, being similar in overall configuration, are not randomly distributed in 
this m«-dimensional image space and therefore they can be described by a relatively low 

10 dimensional subspace. Using principal component analysis (PC A), the vectors that best 
account for the distribution of face images within the entire image space can be found. PCA 
involves determining the principal eigenvectors of the covariance matrix corresponding to 
the original face images. These vectors define the subspace of face images, often referred to 
as the face space. Each vector represents an //i-by-zz image and is a linear combination of the 

15 original face images. Because the vectors are the eigenvectors of the covariance matrix 
corresponding to the original face images, and because they are face-like in appearance, they 
are often referred to as eigenfaces [4]. 

When an unknown image is presented, it is projected into the face space. In this way, 
it is expressed in terms of a weighted sum of eigenfaces. 

20 In the present embodiments, a closely related approach is used, to generate and apply 

so-called "eigenblocks" or eigenvectors relating to blocks of the face image. A grid of 
blocks is applied to the face image (in the training set) or the test window (during the 
detection phase) and an eigenvector-based process, very similar to the eigenface process, is 
applied at each block position. (Or in an alternative embodiment to save on data processing, 

25 the process is applied once to the group of block positions, producing one set of eigenblocks 
for use at any block position). The skilled man will understand that some blocks, such as a 
central block often representing a nose feature of the image, may be more significant in 
deciding whether a face is present. 

30 Calculating Eigenblocks 

The calculation of eigenblocks involves the following steps: 
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(1) . A training set of N T images is used. These are divided into image blocks each of 
size m x n. So, for each block position a set of image blocks, one from that position in each 
image, is obtained: \l 0 ' } /= r , . 

(2) . A normalised training set of blocks {/' JJ , is calculated as follows: 
Each image block, /„' , from the original training set is normalised to have a mean of 

zero and an L2-norm of 1, to produce a respective normalised image block, I' . 
For each image block, I J , t = l ..N T : 

_ I J -mean J J 
-||/;-mean_/;| 

where mean J J = XZ^t'^l 

and \\l 0 ' -mean / 1 = l±±(l 0 l [iJ]-me a n_I 0 'J 

11 ™ \ m j=i 

(i.e. the L2-norm of (/„' - mean_/ c ' )) 
(3). A t rainin g set of vectors {x'}^ is formed by lexicographic reordering of the pixel 
elements of each image block, /' . i.e. Each m-by-n image block, 7' , is reordered into a 
vector, x , of length N=mn. 
15 (4). The set of deviation vectors, D = {x'}%, is calculated. D has N rows and N T 
columns. 

(5). The covariance matrix, S , is calculated: 

S = DD T 

E is a symmetric matrix of size NxAf. 
20 (7). The whole set of eigenvectors, P, and eigenvalues, X t , i = 1,.., N , of the covariance 

matrix, S , are given by solving: 

A = / >r 2P 

Here, A is an Nx N diagonal matrix with the eigenvalues, X, , along its diagonal (in 
order of magnitude) and P is an N x N matrix containing the set of N eigenvectors, each of 
25 length N. This decomposition is also known as a Karhunen-Loeve Transform (KLT). 
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The eigenvectors can be thought of as a set of features that together characterise the 
variation between the blocks of the face images. They form an orthogonal basis by which 
any image block can be represented, i.e. in principle any image can be represented without 
error by a weighted sum of the eigenvectors. 

5 If the number of data points in the image space (the number of training images) is 

less than the dimension of the space (N T <N), then there will only be N T meaningful 
eigenvectors. The remaining eigenvectors will have associated eigenvalues of zero. Hence, 
because typically N T < N , all eigenvalues for which i > N T will be zero. 

Additionally, because the image blocks in the training set are similar in overall 

10 configuration (they are all derived from faces), only some of the remaining eigenvectors will 
characterise very strong differences between the image blocks. These are the eigenvectors 
with the largest associated eigenvalues. The other remaining eigenvectors with smaller 
associated eigenvalues do not characterise such large differences and therefore they are not 
as useful for detecting or distinguishing between faces. 

15 Therefore, in PC A, only the M principal eigenvectors with the largest magnitude 

eigenvalues are considered, where M <N T i.e. a partial KLT is performed. In short, PCA 
extracts a lower-dimensional subspace of the KLT basis corresponding to the largest 
magnitude eigenvalues. 

Because the principal components describe the strongest variations between the face 

20 images, in appearance they may resemble parts of face blocks and are referred to here as 
eigenblocks. However, the term eigenvectors could equally be used. 

Face Detection using Eigenblocks 

The similarity of an unknown image to a face, or its faceness, can be measured by 
25 determining how well the image is represented by the face space. This process is carried out 
on a block-by-block basis, using the same grid of blocks as that used in the training process. 
The first stage of this process involves projecting the image into the face space. 

Projection of an Image into Face Space 
30 Before projecting an image into face space, much the same pre-processing steps are 

performed on the image as were performed on the training set: 
(1). A test image block of size mxnis obtained: J 0 . 
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(2). The original test image block, I B is normalised to have a mean of zero and an L2- 

norm of 1, to produce the normalised test image block, / : 

j _ mean I n 



\I 0 -mean_/ 0 || 



1 J 71 , " ^ 

where mean_/ 0 = ^ I Q [i, j] 



5 and ||J 0 -mean_/ 0 l= ^ J(/ c [i,y]-mean_/J 2 

(i.e. the L2-norm of (/, - mean_/ 0 )) 

(3) . The deviation vectors are calculated by lexicographic reordering of the pixel 
elements of the image. The image is reordered into a deviation vector, x' , of length N=mn. 

After these pre-processing steps, the deviation vector, x, is projected into face space 
10 using the following simple step: 

(4) . The projection into face space involves transforniing the deviation vector, x, into its 
eigenblock components. This involves a simple multiplication by the M principal 
eigenvectors (the eigenblocks), P t ,i = 1,..,M . Each weight y t is obtained as follows: 

y f =P, T x 

15 where P t is the i* eigenvector. 

The weights y,, i = l,..,M, describe the contribution of each eigenblock in 

representing the input face block. 

Blocks of similar appearance will have similar sets of weights while blocks of 
different appearance will have different sets of weights. Therefore, the weights are used here 
20 as feature vectors for classifying face blocks during face detection. 



